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REG-2 (CURVILINEAR REGRESSION) REG-4 (LINEAR REGRESSION) 


DECUS Program Library Write-up DECUS No. 8-321 

REG-2 - CURVILINEAR REGRESSION FOR ANY NUMBER OF DATA 


SUMMARY 

This program computes the values a, b and c for a second degree polynomial equation 
of the form: 

Y = a + bX + cX 2 

where X and Y can represent any number of values. The sum of squares of Y, and the sum 

2 

of the products of X Y and XY are also printed out, which are required to determine the 
sum of squares due to regression. Positive or negative values of X and Y may be used. 

TAPES REQUIRED 

Form of program tape - The program is written in the PDP-8 FORTRAN-D language, 
and is in the source language. 

Form of data tape - The data to be analyzed should be punched onto paper tape, the X 
values first, and then the Y values. Both should be expressed as variables, not integers. 

At the beginning of the program the total number of X or Y variables should be punched, 
in the form of an integer (i.e., if there are 40 pairs of X and Y, punch 40). This value 
represents n in the program tape and enables the program to run for any number of variables. 

OPERATING INSTRUCTIONS 

.FORT 

*OUT-S:REG 2 
* 

* IN—R: 

* t 
T 

*READY 

T 

If the program tape has already been compiled onto the disk, it may be called back into 
core as follows: 


Program tape in high-speed reader 


Data tape in high-speed reader 


. FOSL 

*IN-S:REG 2 
*OPT- 
* T 

*READY 

T 








OUTPUT 


The program prints out the values under: 

A, B, C, PXZY, PXY and PYZ 
STORAGE AND LIMITATION 

X and Y values can be stored with up to 3 digits after the decimal point without losing 
too much precision (error in the range of 1(P) for fairly low numbers of variables, though 
it may be slightly higher for high numbers). 

METHOD 

See "Principles and Procedures of Statistics", p. 338, Steel, R. G. D. and Torrie, J. H. 
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CURVILINEAR REGRESSION 2ND METHOD 
DIMENSION X1(25),Y(25) 

DO 20 1=1,25 
READ 2,101,XI (I) 

FORMAT (/, E) 

CONTINUE 
DO 30 J=l,25 
READ 2,102,Y(J) 

FORMAT (/, E) 

CONTINUE 
READ 2,1 04, K 
FORMAT (/, I) 

SV=0.0 

sxi=0.0 

SX2=0.,0 

SX1Y=0.0 

SY2-=0.0 

SX3=0.0 

SX4=0.0 

SX2Y=0.0 

DO 10 N=1,25 

SY=SY+Y(N) 

SX1=SX1+X1(N) 

SX1 Y=SX1 Y+Xl (N)*Y(N) 

IF(X1 (N) )70,71 ,71 
SX2=SX2+X1(N)**2 
SX3=SX3+X1 (N)**3 
SX4=SX4+X1 (N)**4 
SX2Y=SX2Y+(X1 (N)**2)*Y(N) 

GO TO 100 
X1(N)=-X1(N) 

SX2=SX2+X1 (N)**2 
SX3=SX3-X1 (N)**3 
SX4=SX4+X1 (N)**4 
SX2Y=SX2Y+(X1(N)**2)*Y(N) 

IF(Y(N) )80,81,81 
Y(N)=-Y(N) 

SY2=SY2+Y(N)**2 

CONTINUE 

XlM=SXl/25.0 

X2M=SX2/25.,0 

YhA=SY/25.0 

PX3=SX3-SXl*SX2/25.0 

PX2Y=SX2Y-SX2*SY/25.J0 

PXY=SX1 Y-Xl M*YM*25.J0 

IF(X1M)90,91,91 
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REG 2 cont'd 


X1MP=-X1M 

PX2=SX2-(X1 MP**2)*25.0 
GO TO 200 

PX2=SX2-(X1M**2)*25.0 
!F(SX2)110,111,111 
SX2P=-SX2 

PX4=SX4-(SX2P**2)/25 .0 
GO TO 3000 
PX4=SX4-(SX2**2)/25 
IF(PX3)1?0,121,121 
PX3P--PX3 

D=(PX2*PX4)-PX3P**2 
GO TO 400 
D=(PX2*PX4)-(PX3**2) 

PA=(PXY* PX4- PX2Y* PX3)/ D 
PB=(PX2* PX2Y- PX3* PXY)/(D*]0.0** K) 
C=PB 
B=PA 

A=YM-(C*X2M*10.0**K)-(B*X1 M) 

A=A*10.0**K 

TYPE 103, A, B, C 

FORMAT^, ,, A",E / "B ,, / E / "C",E,) 

TYPE 105, PX2Y, SY2, PXY 

FORMAT (/, "PX2Y", E, M SY2", E, "PXY", E) 

END 
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REG-4 - LINEAR REGRESSION FOR SIMPLE LINEAR, EXPONENTIAL 
AND HYPERBOLIC RELATIONSHIPS 


SUMMARY 



This program 

computes the values of: 


(a) 

a and b in the equation 

Y=a + b X 

(b) 

log^ a and log^ b in the equation 

log Y = log a + X log b 
6 © 0 

(c) 

log^ a and b in the equation 

log Y = log a + b log X 
6 © © 


Any number of pairs of X, the independent variable, and Y, the dependent variable, may 
be used and they may be either positive or negative. Any number of regressions may be 
calculated consecutively, using any or all of the above relationships. 


TAPES REQUIRED 


Form of program tape - The program is written in the PDP-8 FORTRAN-D language, and is 
in the source language. 

Form of data tape - The data tape contains the X, Y values for the first regression 
calculation in the form: 


x, Y , 


X 2 Y 2 


X Y 
n n 


If more than one regression calculation is required, subsequent sets of X and Y follow. 


OPERATING INSTRUCTIONS 





.FORT 

*OUT-S:REG 4 
* 

*IN-R: 

* *|S 

*READY 


Program tape in high-speed reader 


Data tape in high-speed reader 


The program will then pause for entry, on the keyboard, of the sequence: 


p q n "Return" 


where p is the number of regressions to be calculated; 

q is 0 for linear regression with no transformation 
1st regression 1 for log^ transformation of Y (exponential) 

2 for log^ transformation of X and Y (hyperbolic) 

n is number of pairs of X, Y variables 
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if p 1 the program will pause for entry of the relevant q n at the beginning of each 
regression calculation 


OUTPUT 

For each regression the program prints the form of the relationship, X, Y, SSX, SSY, SXY, 
2 

a, b, r, r , regression SS, residual SS, mean square deviation, and number of pairs of 
X, Y variables. 

Where tranformations have been carried out, the appropriate transformed parameters 
will be output. 


METHOD 


Snedecor "Statistical Methods". 
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REG 4 


L 

C SYKES REGRESSION Y ON X WITH LOG TRANSFORM OPTIONS REG4 

N5=J0 

1J00 FORMAT (/, E, E) 

101 FORMAT (I, I) 

ACCEPT 101, N1 

128 XBAR=0. 

SS X=0. 

YBAR=0. 

SSY=0. 

SSXY=0. 

ACCEPT 101, N2, N3 
DO 1 J=1, N3 
EJ=J 

READ 2,1^0,X, Y 
IF (N2-l)l10,105,1,07 
105 Y=LOGF(Y) 

GO TO 110 
1,07 Y=LOGF(Y) 

X=LOGF(X) 

110 SSX=SSX+(X-XBAR)*(X-XBAR)*(1. -1 ./EJ) 

SSY=SSY+(Y-YBAR)*(Y-YBAR)*(1 .-1 ./EJ) 
SSXY=SSXY+(X-XBAR)*(Y-YBAR)*(1 .-1 ./EJ) 

XBAR=XBAR+(X-XBAR)/EJ 
YBAR=YBAR+(Y-YBAR)/ E J 
1 CONTINUE 

IF (N2-1) 129,132,133 

129 FORMAT {/,/, "Y=A+B*X", E) 

TYPE 129 

GO TO 134 

132 FORMAT (/,/, "LN(Y)=LN(A)+LN(B)*X", E) 

TYPE 132 

GO TO 134 

133 FORMAT (/,/, "LN(Y)=LN(A)+B*LN(X)", E) 

TYPE 133 

134 B=SSXY/SSX 
A=YBAR-B*XBAR 
R=SSXY/SQTF(SSX*SSY) 

RSQ=R*R 

SSTO=SSXY*SSXY/SSX 

SSFR=SSY-SSTO 

SDEV=SSFK/(EJ-2.) 

116 FORMAT (/,"XBAR",E) 

117 FORMAT (/, "SUM SQUARES X", E) 

118 FORMAT (/, "SUM SQUARES Y", E) 

119 FORMAT (/, "YBAR", E) 

141 FORMAT (/, "SUM CROSS PROD", E) 
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REG 4 cont'd 


120 FORMAT (/, "A",E) 

121 FORMAT E) 

122 FORMAT {/, "R", E) 

123 FORMAT (/, "RSQ", E) 

124 FORMAT (/, "SS TO REG ", E) 

125 FORMAT (J, "SS FROM REG ", E) 

1 26 FORMAT (/, "MEAN SQ DEV", E) 

127 FORMAT (/, "NO PAIRS OBS", E) 
TYPE 116,XBAR 

TYPE 117, SSX 
TYPE 118, SSY 
TYPE 119, YBAR 
TYPE 141,SSXY 
TYPE 120, A 
TYPE 121, B 
TYPE 122, R 
TYPE 123, RSQ 
TYPE 124, SSTO 
TYPE 125, SSFR 
TYPE 126,SDEV 
TYPE 127, EJ 
N5=N5fl 

IF (N5-N1) 128,2,2 
2 CONTINUE 

END 


8 




